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QED vacuum fluctuations and induced electric dipole moment of the neutron 
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Quantum fluctuations in the QED vacuum generate non-linear effects, such as peculiar induced 
electromagnetic fields. In particular, we show here that an electrically neutral particle, possessing 
a magnetic dipole moment, develops an induced electric dipole- type moment with unusual angular 
dependence, when immersed in a quasistatic, constant external electric field. The calculation of this 
effect is done in the framework of the Euler-Heisenberg effective QED Lagrangian, corresponding to 
the weak field asymptotic expansion of the effective action to one-loop order. It is argued that the 
neutron might be a good candidate to probe this signal of non-linearity in QED. 
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It has been known since the pioneering work of Eulcr 
and Heisenberg that quantum fluctuations in the 

QED vacuum induce non-linearity. This happens in spite 
of QED being a linear theory at the Lagrangian level. 
Given the extreme smallness of these induced non-linear 
effects, this subject has remained in relative obscurity 
for many years. Recent renewed interest in this field 
has been motivated in part by the steady increase in 
the intensity of modern lasers [H with associated peak 
electric fields reaching 10 14 V/m, and envisaged fields of 
10 15 — 10 16 V/m in the near future. Such intense fields 
are approaching the critical value E c ~ 10 18 V/m be- 
yond which unusual properties of the QED vacuum are 
expected to manifest themselves. Among these proper- 
ties one has vacuum birefringence, and non-linear Comp- 
ton scattering j4j. In addition to such time-dependent 
non-linear phenomena, classical charge/current sources 
immersed in external quasistatic fields also generate un- 
expected electromagnetic field configurations induced by 
QED vacuum fluctuations. In fact, we have shown re- 
cently [f| (see also @) that an electrically charged (non- 
magnetic) particle immersed in a constant external mag- 
netic field generates several electric and magnetic multi- 
poles with unusual angular dependences. Among these 
there is an electric multipole term, independent of the 
external magnetic field, which corresponds to the correc- 
tion to the Coulomb field due to quantum fluctuations in 
the QED vacuum. More importantly, though, the leading 
term of the induced magnetic field is due to an induced 
magnetic dipole moment. For accessible values of the 
relevant parameters, this induced effect is comparable to 



the magnetic dipole moment of the nucleon. 

Given the intrinsic relationship between electric and mag- 
netic phenomena in Electrodynamics, the reciprocal ef- 
fect is to be expected, i.e. that a magnetic dipole im- 
mersed in a constant electric field might generate an elec- 
tric dipole moment. In this paper we show that this is in 
fact the case, and argue that the neutron might serve as a 
probe of this induced non-linear effect in QED. While, in 
principle, these non-linear effects could be studied within 
the full QED theory, their complexity calls for a simpli- 
fying approach such as that of a corresponding effective 
theory. For instance, looking at photon-photon scatter- 
ing mediated by fermionic loops, and in the kinematical 
regime where the photon momenta arc much smaller than 
the fermion rest masses, there follows an effective QED 
action. In the case of quasiatatic electromagnetic fields, 
the first non-linear term in the weak field expansion is 
the Euler-Heisenberg Lagrangian 
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where the omitted terms are of higher order in the ex- 
pansion parameter f, which in SI units is given by 
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with a = e 2 / (Aireohc) the electromagnetic fine structure 
constant, and c the speed of light. The invariants T and 
Q are defined as 



T = 



1 



E 



B 2 ) 



I F F^ v 

4 ^ 



(3) 



2 



cE B 



1 



j. I V 



F^ . 



(4) 



with = dpAy-dvAn and = \t^F p<J . Clearly, 
it is possible to use a different, dimensionless expan- 
sion parameter of the QED effective action in terms of 
e.g. a critical electric or magnetic field for the onset 
of non-linearity. This critical field is roughly given by 
that needed to produce an electron-positron pair from 
the vacuum in a length scale of a Compton wavelength, 
i.e. E c = m 2 e c 3 /he ~ 1.3 x 10 18 Volt/m, where m e and e 
are the mass and charge of the electron, respectively. In 
this case, £ = K(m e /E c ) 4 , where k = 2a 2 e 2 c 7 /(45e 4 ?i), 
or using natural units (h = c = 1) one has k = l/(3607r 2 ). 
However, in the sequel we use £ as the expansion param- 
eter. The constitutive equations which follow in the non- 
linear effective theory are very different from the stan- 
dard ones: D = e E + P, and H = — M, where 

' Mo 

P = M = 0, if quantum vacuum fluctuations are ignored. 
In fact, from the full Lagrangian, Eq.(l), it follows that 
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These equations imply the well known non-linear rela- 
tions P = 2£(4.FE + 7c£B) , and M = 2£(-4c 2 FB + 
7c£E) . In terms of the fields D and H the Euler- 
Lagrange equations of motion reduce to the linear 
Maxwell equations with all the effects of non-linearity 
contained in the constitutive equations, Eqs.(5) and (6). 
Therefore, the equations for D and H can be solved as in 
the usual linear theory, i.e. V-D = jo, — <9D/<9i+VxH = 
j, V • B = 0, V x E + dB/dt = 0. The relations between 
the fields D and H and the electromagnetic intensities E 
and B, Eqs.(5) and (6), can be easily inverted to leading 
order in £ (see @). As shown in [j| it is possible to ob- 
tain a general analytical solution to the inhomogeneous 
equations for D(x) and H(x) in terms of the gradient of 
a scalar function 0(x), and the curl of a vector function 
K(x), both of order 0(C), and calculable in terms of the 
external electric/magnetic sources. One first writes 0] 



D(x) =D(x) M +V x K(x) 



H(x) = H(x) 
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where D(x)m and H(x)m are the solutions in the linear 
theory which satisfy V x Dj\/ = 0, and V • Hjf = 0. The 
vectors V x K(x), and V</>(x) can be obtained using the 
remaining homogeneous equations with the result 
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To isolate the induced (non-linear) effects we define new 
electric and magnetic fields with respect to the (linear) 
Maxwell theory as 



E (x) = E(x) - -D M (x) , 



B(x)=B(x)- Mo Hm(x), 



(11) 



(12) 



with V x 8 (x) = 0, and V • B(x) = 0. Using in Eqs.(ll)- 
(12) the expressions for E(x) and B(x) that result from 
inverting Eqs.(5)-(6), together with Eqs.(7)-(10), it fol- 
lows to leading order in £ 
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The two equations above provide general analytical 
expressions to compute the induced fields in terms of the 
Maxwell fields. The latter can be obtained as usual, once 
the external classical charge and/or current distributions 
are specified. 



We now consider as a source an electrically neutral parti- 
cle of radius a, possessing a magnetic dipole moment, in 
the presence of an external, quasistatic, constant electric 
field Eo. The details of how this magnetic dipole moment 
is produced will play no role in the sequel. In fact, the 
final result for the induced electric dipole moment, gener- 
ated by non-linearity in the presence of E , depends only 
on a and m (and obviously on Eo). Hence, for simplic- 
ity, we consider a current density uniformly distributed 
on the surface of a sphere of radius a 
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This current distribution gives rise to a magnetic dipole- 
type field 

i_> Mo f [3 (m ■ e r ) e r - m] 2m rv a 

Bd = ^{ ~3 e(r-a) + — 6(a-r) 

where m is identified with the magnetic dipole moment 
of the source. A similar result is obtained for a uniformly 
magnetized sphere of the same radius a. Since the central 
expressions, Eqs.(13)-(14), were derived assuming E = 
cB < E Cl the following constraint follows 
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where r c stands for a critical (minimum) radial distance 
from the source, e.g. if |m| ~ 10~ 26 Am 2 , then it follows 



3 



that r c > 10 fm. We shall return to this bound later on in electric field is 
the conclusions. With D M = en E , and H M = — B d , 

Mo 



the densities Tm and Qm become 
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Choosing the magnetic dipole moment along the z-axis, 
m = |m|e z , and the external electric field Eo in the x-z- 
planc, forming an angle ip with the z-axis, 



E = E (sin ip e x + cos ip e z ) . 



(20) 



it follows that 



(Eq ■ e r ) = Eq (cos ip cos 9 + sin ip sin 9 cos 0) , (21) 



where the angular function in parenthesis is cos 7, where 
7 is the angle between the external electric field and the 
radial direction. Concentrating on the induced electric 
field, the integrand in Eq.(13) is given by 
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Substituing this expression in Eq.(13), expanding the 
inverse distance in terms of spherical harmonics, using 
their orthogonality properties, and performing the three- 
dimensional integration, the final result for the induced 
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where |x| = r. The first term in brackets above, with 
the l/|x| 2 behaviour, corresponds to the induced elec- 
tric dipole moment. For the case in which the external 
electric field E is parallel to the magnetic dipole mo- 
ment along the z-axis, i.e. for -0 = 0, the induced electric 
dipole field is given by 
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The induced electric dipole moment, located at the ori- 
gin, is then given by 
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When m and E are not parallel, the term proportional 
to sin-0 in Eq.(23) leads to an explicit dependence on the 
azimuthal angle <fi, viz. 
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where the induced electric dipole moment is now 
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which reduces to Eq.(25) when ip = 0. This induced elec- 
tric field is of the electric dipole- type in its radial l/|x| 3 
dependence, but it has a manifestly peculiar angular de- 
pendence. For instance, along the z-axis, and unlike a 
standard electric dipole field, it has a non-zero compo- 
nent along eg that depends on the azimuthal angle (p. 
It also has a non-zero component along the direction of 
e^, as may be appreciated by writing the induced electric 
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field in spherical coordinates, i.e. 
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This feature should help to differentiate this induced ef- 
fect from other predictions of an intrinsic electric dipole 
moment of the neutron in extensions of the Standard 
Model [gj. To arrive at an estimate of the size of the 
induced electric dipole moment of the neutron, we con- 
sider first the case ip = 0, i.e. the external electric field Eo 
parallel to the z-axis, so that the induced electric dipole 
field has the usual angular dependence. Expressing |E | 
in V/m, and the radius a in fermi, Eq.(25) gives 
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A reasonable estimate of the neutron radius can be ob- 
tained from its mean squared radius [9»] in which case 
a ~ 0.35 fm . Current experimental sensitivity [l(]| is at 
the level of 10 _26 |e|cm, which could improve by two or- 
ders of magnitude in the SNS experiment at Los Alamos 
PH . in the CryoEDM experiment at ILL [12], or in the 
nEDM experiment at PSI [l3[. However, current experi- 
mental setups involve an external magnetic field, around 
which the neutron spin precesses, and an external electric 
field to probe the electric dipole moment. This external 
electric field, though, is at present well below the one re- 
quired to bring |p|/at_d to an observable level. In fact, 



one would need fields as high as |Eo| — 10 7 — 10 8 (V/m). 
Such intense fields are indeed present in some crystals at 
the level of |Eo| — 10 10 V/m, but would require a com- 
pletely different setup [14j. In this case, with a ~ lfm, 
the order of magnitude of the induced electric dipole mo- 
ment would be IpI/jvzj — 10~ 23 |e|cm. While this range 
is within current sensitivity, it is crucial to change the 
experimental setup so that the external electric field is 
as high as indicated above. Turning to the more general 
case tj) ^ 0, the order of magnitude of y>\ind does not 
change much, but the peculiar angular dependence of the 
electric field offers some hope of improving the chance of 
observability. To aid in the design of the experiment we 
point out that, from Eq.(27), the induced electric dipole 
moment has a component along the direction of Eo, and 
a component along the x-axis. Hence it will interact with 
Eo and change the Larmor frequency of the neutron mag- 
netic moment around the external magnetic field. 
In connection with the bound, Eq.(17), which follows 
from the weak field approximation, it would be satisfied 
quite easily for the neutron if the experiment is performed 
using a crystal, since in this case r c » 10/m. We men- 
tion in closing that we have computed the induced mag- 
netic field from Eq.(14). The leading term in this field is 
of the magnetic dipole type, and thus it is a correction 
to the field produced by the source (neutron). However, 
its magnitude is many orders of magnitude smaller than 
the latter, so that it can safely be ignored. 
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